Phase Diagram And Adsorption-Desorption Kinetics Of CO On i?w(0001): 
Present Limitations Of A First Principles Approach 
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A lattice gas model is used to study the equilibrium properties and desorption kinetics of CO on 
Ru(0001). With interactions obtained from density functional theory (DFT) the phase diagram and 
temperature programmed desorption (TPD) spectra are calculated up to a coverage of 1/3 ML using 
top sites only. For coverages beyond 1/3 ML hollow sites are included. Good agreement is obtained 
between experiment and theory for coverages below 1/3 ML using top sites only. When including 
hollow sites, DFT calculations fail in predicting the correct binding energy differences between top 
and hollow sites giving disagreement with TPD, low energy electron diffraction (LEED) and heat 
of adsorption experiments. 
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I. INTRODUCTION 

Adsorbate systems with commensurate structures 
can be described successfully with a lattice gas model 
with one or more adsorption sites. The energetics are 
contained in a few parameters such as binding energies, 
vibrational frequencies and lateral interactions. In a 
"complete" theory these parameters would be obtained 
from first principles quantum mechanical calculations, 
for instance based on density functional theory. This has 
been demonstrated successfully for oxygen on Ru(0001) 
for both equilibrium and desorption kinetics P, Q . An 
attempt to repeat such an approach for CO/Pt(lll) 
has been shown to fail 3], most importantly due to 
the fact that DFT does not produce the right energetic 
order of the binding sites for many transition and noble 
metal substrates [H, E3- Nevertheless a satisfactory 
explanation of all available equilibrium and kinetic data 
could be achieved by treating a minimal set of the lattice 
gas parameters as phcnomcnological fitting parameters 
In this paper we will show that for CO/Ru(0001) 
DFT yields parameters that give a good account for 
data below 1/3 ML provided only top sites are included, 
but fails when bridge and hollow sites come into play. 

Let us begin with a quick survey of the experi- 
mental data needed to model this system appropriately; 
details and actual data will be given later. At low 
coverage {9 < 1/3 ML) and temperature (T < 150 
K) CO binds to on-top adsorption sites only resulting 
in a (V3 x ^3)^30° ordered structure at 1/3 ML 
m EH • As the coverage is increased beyond 1/3 



ML three fold hollow sites become occupied forming a 
(2V3 x 2\/3)i?30 o structure at 1/2 ML, with equal pop- 
ulation of top, hep and fee sites, for which two different 
geometries have been proposed 0, 0|. For coverages 
beyond 1/2 ML, a p(7 x 7) structure at 0.55 ML with top 
and bridge sites [ill and a (2^3 x 2V3)i?30° structure 
at 7/12 ML 11. |12J have been suggested. Finally, a 
(5V3 x 5V3)i?30° saturation structure at 49/75 ML 
was observed using He scattering experiments |l3j . A 
summary phase diagram was given by Pfniir et al. |12| . 
Equilibrium isobars and the isosteric heat of adsorption 
|14j have been measured, too, as well as the temperature 
programmed desorption spectra [lH Il5j |. The sticking 
coefficient was measured directly with molecular beam 
scattering as a function of temperature, coverage and 
kinetic energy of the incident beam EH and by coverage 
vs. exposure curves [lfj. Measurements have also 
been done for the CO stretch frequency u sing infrared 
reflection-adsorption spectroscopy (IRAS) j^J 0, Ils| . 
The other frequency modes of the on top species have 
been determined using helium scattering experiments 
[l3l | and IRAS measurements [T^| . 
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In the next section we briefly introduce the lattice gas 
model, the DFT setup as well as the methods used to 
calculate the equilibrium and kinetic data. In section 3 
we consider the modeling of CO/Ru(0001) with top sites 
only up to 1/3 ML and in section 4 we show and discuss 
the modeling of this system with top and hollow binding 
sites up to 1/2 ML. The paper ends with a discussion to 
what extent present DFT methods can be used to model 
CO/Ru(0001). 
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II. THEORETICAL METHODS 
A. Lattice gas formalism 

To set up a lattice gas model for CO/Ru(0001) we 
require a Hamiltonian, which we write down for top sites 
only as: 

H = -(V + k B Tln(q 3 q int ))J2U (1) 

i 

ti+b 



n % a n 



i a, 6 



with obvious generalizations to a multi-site system 0, |j| . 
Here the sum on i exhausts all lattice cells; i + a n , 
i + a and i + b label neighboring cells and we have 
introduced occupation numbers ti — or 1 depending 
on whether a site in cell i is empty or occupied. Vq 
represents the depth of the adsorption potential well 
with reference to the gas phase molecule. Moreover, 
Vi, V2 and V3 denote first, second and third nearest 
neighbor interactions respectively. Vtrio includes Vu, 
Vbt, Vu which denote linear, bent and triangular trio 
interactions involving first and sometimes second nearest 
neighbor interactions Q, Furthermore, 53 = q 2 q xy 
is the vibrational partition function of the adsorbed 
molecule for its center of mass vibrations with respect to 
the surface with q z being the component perpendicular 
to it. Likewise, q xy is the partition function for the 
motion parallel to the surface. We have also allowed 
for the fact that the internal partition for rotations 
and vibrations of an adsorbed molecule is changed from 
its free gas phase value, Z int to gj nt , if some of the in- 
ternal degrees of freedom are frozen out or frustrated [2(| . 

We determined the temperature-coverage phase di- 
agram by calculating the corresponding (-\/3 x a/3)-R30° 
order parameter |2l| . Second order phase transitions 
were defined at a given temperature by the inflection 
point of the order parameter [52, as a function of 
the chemical potential. First order phase transitions 
were marked when it was discontinuous [23L |24| . To 
calculate the resulting TPD spectra we used the theo- 
retical framework described in [2(J • The desorption rate 
depending essentially on the sticking coefficient S(6,T) 
and the chemical potential of the adsorbate fi(8,T) 
(which includes its binding energy and the frequencies of 
the CO molecule in the gas phase and on the surface). 



B. Monte Carlo Methods 

The TPD, isobars and the phase diagram were calcu- 
lated with Monte Carlo methods. Simulations were per- 
formed in both the grand canonical and canonical ensem- 
bles (chemical potential and coverage specified, respec- 
tively) using the Metropolis algorithm (with spin flip and 



infinite Kawasaki dynamics) |25 ; ]. Equilibration times of 
the order 2 14 up to 2 18 Monte Carlo sweeps were allowed 
for each coverage or chemical potential point. When 
dealing with top sites only, we initialized the system in 
a (v3 x \/3)R30° structure for the first temperature- 
coverage, temperature-chemical potential point, respec- 
tively. For coverages beyond 1/3 ML, we occasionally ini- 
tialized the system with a p(2 x 2) structure (with equal 
populations of top and hep sites at 1/2 ML) or with a 
clean substrate. Thereafter, to fasten equilibration, ini- 
tialization of the system occurred with the final config- 
uration of the previously calculated point. When calcu- 
lating the phase diagram, we performed averages over a 
maximum of five independent samples. For the canoni- 
cal ensemble, the particle insertion method of Widom j2|J 
was used in order to calculate the chemical potential with 
Monte Carlo methods; a method which was developed 
to be applicable especially for the interaction parame- 
ters derived from DFT calculations. To allow adsorption 
on hollow and top sites we used three interpenetrating 
42 x 42 or 60 x 60 lattices, with periodic boundary condi- 
tions. This lattice size was chosen to allow commensurate 
low-coverage ordered structures with periods of 3 and 6. 



C. Density Functional Theory setup 

We have calculated all necessary interaction parame- 
ters with the Vienna ab initio simulation package VASP 
[271 l28l l29j . a plane-wave DFT program, which is based 
on the projector augmented wave method |2jj. A cut- 
off energy for the expansion of the plane waves of 400 
eV was found to be sufficient for an accurate descrip- 
tion. For exchange and correlation the generalized gra- 
dient approximation (GGA) according to Perdew et al. 
[3l| was applied. The surface was modeled in p(2 x 2) 
and p(3 x 3) cells at the theoretical lattice constant of 
a=2.725 A and c/a=1.579, with slabs consisting out of 
6 layers, of which the uppermost two layers were struc- 
turally optimized where mentioned in the text. The Bril- 
louin zone was sampled by a grid of (3 x 2 x 1) k-points. 

The adsorption energy for an isolated CO molecule 
(—Vq) is defined as the adsorption energy for 1/9 ML 
coverage (E^ -1 ^ 9 ) derived from the total energies of the 



adsorbate/substrate system (E 1 



CO/Ru^ 
total 



the bare surface 



(*&) and the free CO molecule (EgfJ: 



V - t,„ - hi 



CO/Ru 
total 



pRu 
^ total 



Efofal (2) 



The interaction parameters where determined by calcu- 
lating the energy of the system for ten configurations at 
various coverages between 1/9 and 1 ML [lj . The energy 
of the system at a coverage of 1 ML was calculated twice: 
once with the (3 x 3) unit cell and again with a (2 x 2) 
unit cell. Consequently, all our interaction sets had 
more ordered structures than lattice gas Hamiltonian 
parameters so that we have fitted, using a least squares 
procedure, for the resulting interactions. Maximum 
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deviations from the calculated energies of the ordered 
structures with respect to what was obtained from the 
interaction parameters, were never greater than 5 meV. 



III. MODELING WITH ON-TOP SITES AT 
LOW COVERAGE 

CO/Ru(0001) has already been modeled successfully in 
the past using a lattice gas model with top sites only. The 
interaction and binding energy parameters were obtained 
phenomenologically by matching the phase diagram [33 , 
the TPD spectra and the isosteric heat of adsorption [33[ . 

In this work, three strategies were used to obtain these 
energies: 

• In the first (method A) the top two layers of the 
bare substrate were completely relaxed and then 
kept fixed at these positions for the description of 
CO adsorption. The geometry of the adsorbed CO 
molecule was optimized only for the 1 /9 ML config- 
uration. For higher coverages, the same CO bond 
length and adsorption height was used. 

• In method B the position of the CO molecules as 
well as the first two layers of the substrate were 
allowed to relax for all of the ordered structures. 

• For comparison, we have also looked at the result- 
ing interaction parameters when fixing laterally the 
CO molecules for the structure at 2/9 ML (which 
is the only configuration for which lateral forces on 
the CO molecule are symmetrically allowed) with 
all the other ordered structures fully relaxed (the 
CO molecules and the first two layers of the sub- 
strate). This third approach (method C) was used 
for O/Ru(0001) 0. 

In Table 1, we show the resulting energies using all 
of the above three methods together with the standard 
deviation a describing the average deviation between the 
calculated energies and the parametrization. We used 
the minimum number of parameters so as to have an 
acceptable value of a. These interaction sets show that 
the second neighbor interaction between top sites, V2, is 
repulsive for the relaxed calculations. This is clearly at 
odds with the experimental phase diagram which implies 
a coexistence between the (-\/3 x v / 3)-R30° structure and 
a lattice gas, and requires V2 attractive. In addition, the 
calculated TPD spectra using transfer matrix techniques 
with method B and C give low initial-coverage spectra 
peak positions that do not correspond to experiment and 
are much too broad. On a more positive note, the first 
neighbor interaction between ontop sites for method B 
(cf. Table P) is closer to what was obtained by previous 
phenomenological models of this system using only one 
binding site [33j. 



Interestingly, for method C V2 becomes negative (—6 
meV) if the same set of parameters (Vi, V2, Vh, V u ) is 
used as for method A. However, in this case the standard 
deviation increases significantly to 19 meV. 

Because of this obvious discrepancy for the relaxed cal- 
culations, we restrict the further analysis, to the results 
obtained within method A. The reason for this behav- 
ior can be found probably in local relaxations for specific 
configurations used for the determination of the param- 
eters, which are not representative for the interaction of 
CO on Ru(0001) in general. For the unrelaxed setup in 
method A, the CO molecule and the surface is treated 
completely identically in all configurations, entirely in 
the spirit of the lattice gas like interaction. For the re- 
laxed setups, however, the number of input configura- 
tions would probably have to be increased such that re- 
laxation effects, which are only characteristic of one spe- 
cific arrangements of molecules, but not of the interaction 
between molecules itself, are averaged out. 

The resulting phase diagram and TPD spectra with 
the parameters of method A are given in Fig^and show 
reasonable agreement with experiment provided we lower 
the binding energy of a single CO molecule by -174 meV 
(our unadjusted binding energy places the model curves 
55 K too high). The experimental sticking coefficient was 
used 0| . Even with this adjustment there is still some 
disagreement for the spectra with an initial coverage of 
1/3 ML. One way to correct for this would be to allow 
for hollow sites in our model for which there would be a 
small but finite number around these temperatures and 
coverages. 

The phase diagram calculation does not depend on the 
binding energy of the top sites. It is dominated by two 
phases, the (\/3 x \/3)R30° structure, which saturates at 
1/3 ML and a non-ordered lattice gas (l.g.) phase. To 
access the high temperature part of the phase diagram 
experimentally would require pressures up to mbar, 
which is not accessible with LEED It has in fact been 
predicted that T c must be at least 550 K at 1/3 ML. This 
is in agreement with our calculations, which predict an 
order disorder temperature of 1000 K at 1/3 ML under 
the restriction of on top adsorption only. We remark 
that this value will be considerably lowered if there is 
a spillover into other types of binding sites for cover- 
ages beyond 1/3 ML, resulting in a continuous reduction 
in temperature (and chemical potential) of the hight of 
the order phase boundary around 1/3 ML. Moreover, the 
value of V2 can be directly read from the experimental 
phase diagram |32| with the value of the tricritical point 
given by « I.O 6V9 . Experimentally, the tricritical point is 
around 150 K |l2j| for which there is some evidence from 
slow diffusion of the CO molecules 8] . This could be the 
source of the small discrepancy between our calculations 
and experiment for the first order transition points. 

From these comparisons with experiment we can 
conclude that DFT does reasonably well in at low 
coverage (9 < 1/3 ML) and since even in this coverage 
regime hollow sites come into play at high temperatures 
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TABLE I: Interactions energies (eV) derived within the DFT 
framework with three different strategies (A,B,C, cf. text) 
for the CO/Ru(0001) system for top sites together with the 
corresponding standard deviation (a). E top denotes the bind- 
ing energy of the top site, the definition of the interaction 
parameters is given in the text. 



Method 


Etop 


Vi 


v 2 


v 3 


Vu 


V bt 


Vu 


a (meV) 


A 


-1.796 


0.230 


-0.010 






0.010 


-0.048 


3 


B 


-1.948 


0.137 


0.002 


0.022 


0.104 


0.066 


-0.231 


2 


C 


-1.948 


0.263 


0.002 


0.022 


-0.021 


0.004 


-0.042 


2 



(T > 150 K) it does even better at low tempera- 
ture. Similarly, the adsorption isobars at 1.3xl0 -4 mbar 
and 1.3xl(T 6 mbar are limited to 1/3 ML for T > 400 K. 

Moreover, one may question the relevance of incor- 
porating trio interactions for coverages below 1/3 ML. 
Indeed, with such a large nearest neighbor interaction 
these trios will practically never occur below 1/3 ML, 
even at desorption temperature and despite the attrac- 
tive triangular trio interaction. However, fitting the 
calculated energies with V\ and V% alone using method 
A results in V x = 0.216 eV and V 2 = -0.060 eV with 
(7 — 4 meV. This would give a phase diagram similar 
to that given in Fig. ^ but with a tricritical point that 
one can estimate [32| to be 75 K instead of ~ 120 K as 
shown in Fig. ^ For coverages beyond 1/3 ML, it is 
quite plausible that trio interactions become important. 
Indeed, this would help to explain the asymmetry in 
the experimental phase diagram. However, as been 
noted before, the asymmetry of the phase diagram can 
be explained either with a spillover into other binding 
sites |22| or with trio interactions 2j, but it is quite 
possible that there is a combination of both. Since from 
the experimental data we know that a spillover must 
occur for coverages beyond 1/3 ML one cannot argue 
from symmetry arguments alone that one must have 
trio interactions. This makes it necessary to analyze 
the experimental data for coverages beyond 1/3 ML 
to consider the relevance of trio interactions for this 
system. With all of these facts in mind, we now proceed 
to model the adsorption behavior at higher coverages. 



IV. MODELING WITH ON-TOP, HCP AND 
FCC SITES 

The interaction energies between sites of the same type 
for fee and hep sites were calculated with the same pro- 
cedure as for the top sites. For the same reason as for 
the low coverage regime, we present here the results of 
method A only. In addition, we show in Table |H] the 
adsorption energies per CO molecule for all investigated 
structures. These results are in line with previous DFT 
calculations that predict that the top sites are the most 
favorable binding site at low coverage [34l l35l | . We note 
that the energy for the (1 x l)-CO structure using a (2 x 2) 
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FIG. 1: (a) Phase Diagram up to 1/3 ML in the temperature- 
coverage plane using the lattice gas model with top sites only 
on a 60x60 lattice (interactions from method A, cf. Table 1). 
The crosses are the experimental points fl^ . the circles indi- 
cate second order phase transitions, the triangles first order 
transitions between lattice gas (l.g.) and (\/3x \/3)-R30 o (v3) 
phases, (b) The corresponding phase diagram in the chemical 
potential-temperature plane, (c) The predicted temperature 
programmed desorption spectra using the experimental stick- 
ing coefficient [lf^ and a heating rate of 5 K/s at initial cov- 
erages of 0.33, 0.25 and 0.15 ML with the top site only model 
(solid), the experimental spectra at 0.50, 0.33 and 0.22 ML 
(dashed) [Tsl ] and the spectra with the top and hollow site 
model at 0.50 and 0.33 ML (dotted). The theoretical curves 
have had their binding energy adjusted by -174 meV in order 
to aid the comparison with experiment. 
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TABLE II: Calculated adsorption energies (eV) for CO/Ru(0001) per CO molecule with method A (cf. text). 



site 














P «=3/4 


^0,(2x2) 


p9=l 

^,(3x3) 


top 


-1.796 


-1.686 


-1.793 


-1.827 


-1.572 


-1.452 


-1.371 


-1.172 


-1.176 


hep 


-1.750 


-1.635 


-1.774 


-1.755 


-1.509 


-1.389 


-1.281 


-1.118 


-1.127 


fee 


-1.672 


-1.557 


-1.660 


-1.676 


-1.426 


-1.323 


-1.231 


-1.097 


-1.104 



and a (3 x 3) unit cell differ by 7 meV for fee sites and 9 
meV for hep sites. When determining the corresponding 
interactions, we needed to include linear trios as well as 
third nearest neighbor interactions to obtain a reasonable 
value of a. Thus, we have also included them for the top 
sites for consistency. The resulting interaction parame- 
ters are shown in Table ITTT1 For each interaction between 
different types of sites we used the calculated energy per 
CO molecule (-^structure) for a single co-adsorption config- 
uration in the p(3 x 3) cell. Each interaction was deduced 
by: ^structure = [E a + E b + V a b)/2, where E a and E b 
refer to the adsorption energy of the two different bind- 
ing sites (Table |nj and V a b is the parameter describing 
the interaction between the two sites of different type. 
The calculated frequencies at low coverage (1/9 ML) are 
compiled in Table ITVl The top site frequencies agree rea- 
sonably well with experimental data and previous DFT 
calculations for 1/3 ML Until very recently [Isj . 

only one peak in the typical on-top region was observed 
by IRAS experiments for coverages exceeding 1/3 

ML, indicating either strong coupling and/or a very low 
infra red intensity of the molecules adsorbed in non atop 
sites. On the other hand, hollow and bridge sites have 
been inferred from Helium scattering experiments, but 
the precise values of their frequencies cannot be deter- 
mined at the present time |l3j . However, our calculated 
stretch frequencies for hep and bridge sites are in good 
qualitative agreement with previous DFT calculations at 
1/3 ML that have their stretch frequencies at 1800 cm -1 
and 1885 cm" 1 respectively |34| . 

As can be seen from Table IIIII one has such a small 
binding difference between top and hep sites that at des- 
orption temperatures (T > 300 K) a significant (greater 
than 0.1 ML) population of hep sites below 1/3 ML oc- 
curs, completely at odds with experiment. The resulting 
isobars from these interaction parameters and frequen- 
cies |33| is shown in Figure 2a, with an adjustment of 
the binding energy of -174 meV (the same shift needed 
for the TPD spectra). At first glance, the curves seem 
to agree well with the experimental data points included 
in the figure. However, a closer investigation reveals, 
that the inflection at 1/3 ML is not captured correctly in 
the calculation. This disagreement is more clearly seen 
when the resulting isosteric heat of adsorption is calcu- 
lated from these isobars using the ASTEK software [3(| : 
the calculated heat of adsorption does not rise at 1/3 
ML unlike experiment which rises from ~ 155 kJ/mol to 
over 180 kJ/mol at 1/3 ML as seen in Figure 2b. This is 
a consequence of having a significant population of hep 
sites below 1/3 ML. It has been argued [33 that the ex- 



perimental value of the heat of adsorption is too high by 
about 8 kJ/mol at 1/3 ML, but this does not improve 
significantly the agreement with our calculations. A sim- 
ilar disagreement occurs for the TPD spectra (Figure^) 
when using the same sticking coefficient as with top sites 
only: The absence of a significant rise in the heat of ad- 
sorption at 1/3 ML is reflected in a shallower minimum 
in the desorption rate around 440 K for initial coverages 
greater than 1/3 ML. Correspondingly, the calculated 
rate, at around 450 K, does not peak as sharply as in ex- 
periment. (The small spike in the desorption rate at 430 
K results from a small change in the chemical potential 
at this coverage (0.36 ML). An examination of the cal- 
culated sticking coefficient here shows that this spike 
will be eliminated in a full model calculation.) Moreover, 
the binding difference between fee and top sites is much 
larger than between hep and top sites. This excludes a 
(2-\/3 x 2\/3)i?30 o structure, for which equal population 
of hep and fee sites has to occur at 1/2 ML. 

We also remark that the inflexion point at 1/2 ML for 
both pressure isobars is due to a (4x2) structure with 
equal populations of hep and top sites. This ordered 
structure is reflected in the rise of the heat of adsorp- 
tion and its subsequent sharp drop to about 50 kJ/mol 
around 0.55 ML. In the TPD spectra the peak at 390 
K confirms the presence of the (4x2) structure. How- 
ever, because of fluctuations in the calculated chemical 
potential at this coverage and temperature we have esti- 
mated the height of the peak to have an uncertainty of 
20 %. On a more positive note, DFT calculations seem 
to agree with LEED experiments that suggest that the 
hollow sites have a binding energy between that of the 
top sites and the bridge sites, the latter been calculated 
to have a binding difference with respect to the top sites 
of 147 meV (obtained from the adsorption energy per 
CO molecule at 1/9 ML giving —1.649 eV). Finally, we 
note that the third nearest neighbor interaction between 
top and hollow sites is attractive which does favor the 
(2y/3 x 2\/3)-R30 o structure, but which is not sufficiently 
large to give it. 



V. CONCLUSIONS 

The phase diagram and TPD spectra were calculated 
from first principles and, apart from a difference in over- 
all binding energy, agree well with experimental findings 
up to 1 /3 ML if only top sites are included in a lattice 
gas model. To model the system beyond 1/3 ML we 
included top and three fold hollow sites. The resulting 
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TABLE III: On-site (Vb) and interactions energies (Vi,V2,V3,Vit,Vbt,Vtt) in (eV) to model the CO/Ru(0001) system with top, 
hep and fee sites using method A together with the corresponding standard deviation (a). In the second part interaction 
energies between sites of different type are compiled. 



site 


Vo 


site 


Vi 


v 2 


Vz 


Vu 


V bt 


V u 


a (meV) 


top 


1.796 




0.220 


-0.010 


0.001 


0.008 


0.014 


-0.061 


2 


hep 


1.750 


0.046 


0.240 


-0.001 


-0.008 


0.034 


0.006 


-0.102 


3 


fee 


1.672 


0.124 


0.225 


-0.001 


0.004 


0.015 


0.003 


-0.088 


3 



sites 



distance 
top-hep 
top-fee 
hep-fee 



Vi 
73/3 



v 2 

2V3/3 
0.049 
0.046 
0.044 



V 3 



V21/3 
-0.017 
-0.014 
0.004 



TABLE IV: Frequencies (cm" 1 ) at 1/9 ML coverage obtained 
within the DFT framework (method A) together with exper- 
imental results from Refs. [T3. Hg|. 



Site 




Vy 


v z 




Vrot 


top 


46.0 


40.4 


392.7 


1989 


378.7 


exp. (1/3 ML) 


46.0 


46.0 


445.0 


2025 


413.0 


bridge 


163.0 


41.8 


344.3 


1789 


170.9 


fee 


95.3 


62.9 


300.6 


1759 


158.3 


hep 


149.8 


140.7 


313.3 


1720 


196.1 



binding difference between hollow and top sites was de- 
termined to be too small to reproduce the isobars and 
heat of adsorption and, because of the non-zero binding 
energy difference between hollow sites, does not repro- 
duce the (2-^/3 x 2 v / 3)i?30° as observed in LEED ex- 
periments. These results seem to indicate that present 
GGA functionals are able to reproduce the interaction 
between same adsorption sites, but in order to describe 
the energy differences between different adsorption sites 



further improvement is necessary. This shortcoming is 
not only limited to the CO/Ru(0001) system, but is well 
known for CO adsorption on transition metals in general 
0, 0| . The reason for this is an insufficient description of 
the energy difference between the highest occupied (2-7T*) 
and the lowest occupied molecular orbital (5cx) of the 
CO molecule, which contribute with different intensity 
at various adsorption sites. Hence, inter-site energy dif- 
ferences are affected. Recently, first attempts have been 
undertaken to cure this problem and the results look very 
promising jfj. 
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